########
setwd("D:/入侵/数据/分析2021/4.5 cor")
library(ggplot2)
library(GGally)
theme_customize2=theme_bw()+
  theme(
    panel.grid.major = element_blank(),
    panel.grid.minor = element_blank(),
    axis.line = element_line(colour = "black"),
    legend.title = element_blank())+
  theme(axis.text.x=element_text(face="bold",size=14,angle=0,color="Black"),
        axis.text.y=element_text(face="bold",size=14,angle=0,color="Black"))+
  theme(axis.title.x = element_text(size = 15, color = "Black"),
        axis.title.y = element_text(size = 15,  color = "Black"))

cor_data=read.csv("cor.csv",header=TRUE)

library(Hmisc)
library(corrplot)

###plant
g3=ggpairs(cor_data,
           columns = c(3:7),
           lower=list(continuous="smooth"),
           upper=list(continuous=wrap("cor",size=8)),
           diag=list(dicrete="barDiag"))
g3+theme_customize2  #bw 8*8

###soil

g3=ggpairs(cor_data,
           columns = c(8:13),
           lower=list(continuous="smooth"),
           upper=list(continuous=wrap("cor",size=8)),
           diag=list(dicrete="barDiag"))
g3+theme_customize2  #bw 10*10